DRAFT January 27, 2011 

Preprint typeset using I^T^X style cmulatcapj v. 08/13/06 



EVIDENCE FOR AN FU ORIONIS-LIKE OUTBURST FROM A CLASSICAL T TAURI STAR 

Adam A. Miller 1 , Lynne A. Hillenbrand 2 , Kevin R. Covey 3,4,5 , Dovi Poznanski 1,6,7 , Jeffrey M. Silverman 1 , Io K. 

W. Kleiser 1 , Barbara Rojas-Ayala 3 , Philip S. Muirhead 3 , S. Bradley Cenko 1 , Joshua S. Bloom 1 , Mansi M. 
Kasliwal 2 , Alexei V. Filippenko 1 , Nicholas M. Law 8 , Eran O. Ofek 2- ', Richard G. Dekany 9 , Gustavo Rahmer 9 , 
David Hale 9 , Roger Smith 9 , Robert M. Quimby 2 , Peter Nugent 6 , Janet Jacobsen 6 , Jeff Zolkower 9 , Viswa 
Velur 9 , Richard Walters 9 , John Henning 9 , Khanh Bui 9 , Dan McKenna 9 , Shrinivas R. Kulkarni 2 , and 

Christopher R. Klein 1 

DRAFT January 27, 2011 

ABSTRACT 

We present pre- and post-outburst observations of the new FU Orionis-like young stellar object 
PTF lOqpf (also known as LkHa 188-G4 and HBC 722). Prior to this outburst, LkHa 188-G4 was 
classified as a classical T Tauri star on the basis of its optical emission-line spectrum superposed on 
a K8-type photosphere, and its photometric variability. The mid-infrared spectral index of LkHa 
188-G4 indicates a Class Il-type object. LkHa 188-G4 exhibited a steady rise by ~1 mag over ~11 
months starting in Aug. 2009, before a subsequent more abrupt rise of > 3 mag on a time scale of 
~2 months. Observations taken during the eruption exhibit the defining characteristics of FU Orionis 
variables: (i) an increase in brightness by >4 mag, (ii) a bright optical/near-infrared reflection nebula 
appeared, (iii) optical spectra are consistent with a G supergiant and dominated by absorption lines, 
the only exception being Ha which is characterized by a P Cygni profile, (iv) near-infrared spectra 
resemble those of late K-M giants/supergiants with enhanced absorption seen in the molecular bands 
of CO and H2O, and (v) outflow signatures in H and He are seen in the form of blueshifted absorption 
profiles. LkHa 188-G4 is the first member of the FU Orionis-like class with a well-sampled optical to 
mid-infrared spectral energy distribution in the pre-outburst phase. The association of the PTF lOqpf 
outburst with the previously identified classical T Tauri star LkHa 188-G4 (HBC 722) provides strong 
evidence that FU Orionis-like eruptions represent periods of enhanced disk accretion and outflow, 
likely triggered by instabilities in the disk. The early identification of PTF lOqpf as an FU Orionis- 
like variable will enable detailed photometric and spectroscopic observations during its post-outburst 
evolution for comparison with other known outbursting objects. 

Subject headings: stars: formation - stars: pre-main sequence - stars: individual (LkHa 188-G4, HBC 
722) - stars: variables: T Tauri, Herbig Ae/Be - stars: winds, outflows 



1. INTRODUCTION 

FU Orionis variables are young stellar objects 
(YSOs) that exhibit large-amplitude optical outbursts 
(Am visua l > 4 mag) and remain bright for several 
decades (Herbig 1977). Models in which the radiation 
is dominated by a rapid increase in the disk accretion 
rate have successfully reproduced many of t he observa- 
tional properties of FU Orionis outbursts; see Hartmann 
fc Kenyonj ( |1996[ ) for a review. An alternative interpreta- 
tion invol ving a rapidly rota ting low-gravity star is pre- 
sen ted by Herbig ( 1989[ ) and Herbig, Petrov, & Duemm- 
ler ( 2003 ) . In the first scenario, the cause of the instabil- 
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ity remains uncer tain; possible mechanisms in clude ther- 
mal instabilities (Hartmann & Kenyon 1996) a nd inter- 
actions with com panion stars on eccentric orbits ( Bonnell 
&: Bastien|1992 ). Accretion instabilities are also believed 
to give rise to less dramatic YSO flares. 10 These events 



(e.g. EX Lupi, V1647 Ori, VI 118 Ori Z CMa; see Herbig 
1977] |Herbigl[2007] jAspin et al.|[20T0l |Aspin fc Reipurth 
^0091 TPedele et al.||2007[ |Rettig et a .||2005[ |Lorenzetti 



et~aTl|2007[ |Audard et ai.fl2010j |Szeifert et al.||20101 and 

references therein) are heterogeneous, and it remains to 
be de termined if th ey should be referred to as a single 
class (iHerbigl 120071). These YSO outbursts are never- 
theless aisTmguIsried from FU Orionis-like outbursts in 
that their optical spectra are emission-line dominated 
(whereas FU Ori stars typically show only absorption 
lines during outburst) and the outbursts have shorter 
lifetimes, T out burst ^ 2 yr, th an FU Ori stars which re- 
main in ourburst for decades ( Herbig|[l977 Fedele et al 
20071. 

YSO outbursts are of great interest given their poten- 
tial importance to the broader star and planet formation 

10 In the literature, these events are frequently referred to as 
"EXors", named for the outbursting star EX Lupi. This nomen- 
clature is confusing, however, as EX Lupi is not part of the Orion 
constellation and m any of the so c alled EXo rs seem to have distinct 
prope rties (see e.g., |Fedele et"aL]|2007| |Herbig|[2Tju7| |Covey et "aT] 
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process. Low-mass stars could accrete as much as half 



their final mass d uring FU Orionis outbursts ( Hartmann 



fe Kenyon||1996 l, and the strong winds and outhows 
they are thought to drive likely have a signi ficant im- 



pact on the surrounding in terstellar medium (Croswell, 
Hartmann, fe Avrett||1987[ ). 



One of the major challenges in studying FU Orionis 
stars is the relatively small sample of known examples. 
Although ~20 FU Orionis candidates have been iden- 



tified by Reipurth & Aspin (2010) based on a mix of 
spectroscopic evidence, only a handful (~6-8) of those 
stars have been observed to rise from their pre-outburst 
state to their eruptive state. Only one example, V1057 
Cyg, has had a pre-outburst spectrum taken; it sh owed 
B aimer, Ca II, Fe I, and Fe II lines in emission (Her- 
big 1977). Herbig noted that these features resembled 
those of a T Tauri star, but a lack of detected absorp- 
tion lines precluded a spectral-type classification beyond 
a late-type star based on the red continuum. 

Recently, a new FU Orionis-like variable, associ- 
ated with LkHa 188-G4, was discovered at aj2ooo = 
20 h 58 m 17?00, <5j2ooo = +43°53'42"9 in the direction of 
NGC 7000/IC 5070, also known as the North Amer- 
ica/Pelican Nebula. The source was announced as a 
likely FU Ori-like outburs t on 2010 Aug. 17 UT 11 by 



Semkov & Peneva (2010a), and was later spectroscopi- 
cally c onfirmed as an FU Ori candidate by |Munari et al.| 
(2010). Further photometric observations showed that 
the source became bluer over the course of the outburst 
and that a reflection nebula was e merging around LkHa 
188-G4 ( |Semkov fe Peneva||2010b| . Near-infrared (NIR) 
photometric observation s revealed a ^3 m ag increase in 
the J, H, and K s bands ( Leoni et al.|2010 ). The variable 
was independently discovered by our collaboration dur- 
ing the course of regular monitoring of the North America 
Nebula with the Palomar 48-in telescope. These images 
were automatically reduced, and new variable and tran- 
sient sources were found using the soft ware developed for 



the Palomar Tra nsient Factory (PTF; |Law et al.||2009 
Rau et al.||2009l Bloom et al. 2011, in prep.). On 2010 
Aug. 5 the outburst from LkHa 188-G4 was identified 
and given the name PTF lOqpf. 

Here, we present pre- and post-eruption observations 
of LkHa 188-G4, an FU Orionis-like star, which we argue 
was a classical T Tauri star prior to eruption. Shortly be- 
fore the submission of this article, an independent analy- 
sis^ LkHa 188-G4 was posted on the arXiv by |Semkov| 



et al. ( |2010[ ), in which they also argue that the 2010 erup- 
tion from LkHa 188-G4 is an FU Orionis outburst. In §[2| 
we discuss archival observations of LkHa 188-G4 prior to 
the 2010 eruption. Our new 2009 and 2010 observations 
are presented in § [3j and § [4] contains our analysis of 
those observations. We discuss our conclusions in § [5j 
Throughout this paper we designate the 2010 outburst 
as PTF lOqpf, while we refer to the pre-outburst star as 
LkHa 188-G4. 

2. PRE-OUTBURST OBSERVATIONS OF LkHa 188-G4 

PTF lOqpf is associated with the previously kn own op- 
tical source LkHa 18 8-G4 ([Cohen fe Kuhil|1979|), which 
was later cataloged b y [Herbig fe Bell| ( |1988[ j as HBC 722. 
Cohen & Kuhi ( 1979 ) identified the star as an my = 18.9 



mag emission-line object with line equivalent widths of 
EW Ha = -100.5 A and EW H/3 = -25.8 A; [O I] A6300 
was also present, and NIR photometry was reported. We 
are unable to identify the strong em ission feature near 
4640 A in the |Cohen fe Kuhi| ( |l979| ) spectrum, and the 
authors themselves do not measure or comment on the 
feature. It could b e related to multiply ionized C or N, 
or He II (see, e.g., Sargent & Filippenko 1991), though 
the absence of a known exact wavelength for the feature 
precludes a secure identification. The original spectrum 
is no longer available (M. Cohen, 2010, private communi- 
cation) and we cauti on that the feature may be an error 
in the original data. Cohen & Kuhi (1979) determined 
a spectral type of K7-M0; they also derived an extinc- 
tion estimate of Ay = 3.4 ± 1.2 mag and a bolometric 
luminosity of log(L/L Q ) = 0.43. 

The source is part of a linear c hain of young stellar o b- 
jects that was first designated by Cohen & Kuhi ( 1979 1 as 
the NGC 7000 / IC 5070 LkHa 188 g roup. SyitzerTKAC 
and MIPS observations presented by |Guieu et al. (2009) 
and Rebull et al. (2011) later revealed that the optically 
visible stars are only the surface population of a much 
richer embedded cluster, known as the "Gulf of Mexico" 
cluster due to its location relative to the North Amer- 
ica Nebula. The Tw o Micron All Sky Survey (2MASS; 
Skrutskie et al |2006t 
K s bands on 2000 June 
obtained with the Kitt Peak National Observato ry 0.9 
m telescope was presented by Guieu et al. (2009), while 
broadband ri and narrowband Ha imaging of LkHa 188- 
G4, obtained between 2003 and 2005, is presented in the 
catalog of the INT/WFC Photo metric Ha Survey of the 
North ern Galactic Plane (IPHAS; |Gonzalez-Solares et al. 
2008). We summarize the IPHAS measurements in Ta- 

Owing to crowding, there are no IRAS measurements 
of LkHa 188-G4. Additional IR sources ar e cataloged 
by th e Midcourse S pace Experiment (MSX ; Mill et al. 
19941 and AKARI (Murakami et all 120071). The MSX 



observed the field in the J, H, and 
10. BVI optical photometry 



UT dates are used throughout unless otherwise noted. 



source has a position (aj 2 ooo = 20 58 m 17!59, <5j2ooo = 
43°53'36'/24) which is ~9.2" from the optically derived 
position of LkHa 188-G4. This position has no clear 
counterpart, to within ~3", in the 2MASS catalog, and 
thus the MSX source may represent a blend of emission 
from several stars in the LkHa 188 gro up. There are 
two s ources in the AKARI MIR catalog (Ishihara et al. 
|2010[ ) within 10" of the optical position of LkHa 188- 
G4. Both of these sources are >4.5" from the optical 
position of LkHa 188-G4, however, so we do not include 
them here in our analysis of pre-outburst emission from 
LkHa 188-G4; the photometry of these sources likely 
does not correspond to LkHa 188-G4 alone. AKARI also 
catalogs a source in the far-infrared (FIR) at aj2ooo = 
20 h 58 m 16^809, 6 J200 o = +43°53'41'.'65, which is ~2.4" 
from the optical position of LkHa 188-G4. The closest 
2MASS counterpart to the AKARI FIR source is LkHa 
188-G4, though we cannot rule out the possibility that 
the AKARI flux measurements include a contribution 
from other members of the LkHa 188 group. 

The pre-outburst spectral energy distribution (SED) 
of LkHa 188-G4 is shown in Figure [TJ The pre-outburst 
SED of LkHa 188-G4 exhibits a smooth trend with the 
exception of the two points at ^0.75 and 0.9 /im. These 
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Fig. 1. — Spectral energy distribution of LkHa 188-G4 both be- 
fore (black circles) and during (red squares) the 2010 outburst. Pre- 
outburst NIR photometry is from 2M ASS, while non-simu lt aneous 
optica l and MIR photometry is from IGuieu et a ll (12009k , IRebullI 
|et al.| |2011| l, and |Gon zalcz-Sol ares et al.| |2"oU8| l. F1K photome- 
try from AKARI is shown as open circles because these detections 
may include light from other members of the LkHa 188 group (see 
text and figure |2] l. The solid line shows the model spectrum, from 
Kurucz (1993), of a K7 V star reddened with Ay = 3.1 mag, and 
demonstrates the infrared excess of LkHa 188-G4. The spectral 
slope in the MIR is consistent with the source being a Class II 
YSO prior to outburst. NIR outburst photometry was obtained on 
2010 Oct. 10 with PAIRITEL, while riz photometry was taken on 
2010 Oct. 14 with P60. 

points, corresponding to the IPHAS i band and KPNO / 
band, respectively, are not coeval and differ by a factor of 
^30% in flux, which is similar to the ~0.3 mag variability 
seen in the multi-epoch IPHAS data. 

Direct integration of the pre-out burst SED, excluding 
the AKARI data., results in a total luminosity of ~0.7 L , 
assuming a distance of 600 pc to NGC 7000. A d optin g 
the classification scheme from Lada & Wilking (1984), 
we fit for the MIR spectral index, defined as 

d log XF \ 

a log A 

where a is the spectral index and F\ is the flux den- 
sity at wavelength A. Following a least-squares fit to 
the Spitzer/IRAC data from 3.6 to 8 ^m, we measure 
a = —0.40, indicating that LkHa 188-G4 was a Class 
II YSO prior to its 2010 outburst. Thus, there was no 
evidence of either a dense core of infalling material or sig- 
nificant envelope material in the vicinity of the central 
star. Including the Spitzer/MIPS 24 /im detection re- 
sults in a — —0.77, securely establishing this source as a 
Class II YSO. To demonstrate the MIR excess relative to 
a stellar photosphere, we also show in Figure[T]the model 
spectrum of a solar metallicity K7 V star, reddened with 
Ay = 3.1 mag, which is c onsistent with the inferences 
from Coh en fc Kuhi| ([197 91) . The model s pectrum , taken 

update to the [Kurucz 



from the 



" |Rurucz| \VMS[ 
models, has been normalized to match the flux of LkHa 
188-G4 in the H band. 

Class II YSOs are most co mmonly asso ciated with clas- 
sical T Tauri stars (CTTSs; [Lada||l987[ ). In addition to 
being a Class II YSO, we also know from the multi-epoch 
IPHAS data that LkHa 188-G4 was optically variable be- 
tween 2003 and 2005, years before the observed PTF out- 
burst, with an amplitude of ^0.2-0.3 mag. While some 
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Fig. 2. — P48 ij-band images of PTF lOqpf showing the emer- 
gence of the optical reflection nebula. Each image is centered on 
PTF lOqpf and is 75" on a side. As with many other FU Orionis- 
like variables, the outburst of PTF lOqpf has created an asym- 
metric reflection nebula. Notice that PTF lOqpf is not located in 
the center of the reflection nebula which extends to the south west 
of the pre-outburst point source. All images are registered to the 
same frame with north up and east to the left. All dates are UT. 



CTTSs show variable amplitudes as large as ~2 mag, 
the majority exhibit low-amplitude variability consistent 
with what we observ e in LkHa 188-G4 prior to eruption 
(iGrankm et al.|2007|. The emission-line spectrum, Class 



11 MIR spectral index, and optical variability all point to 
LkHa 188-G4 being a fairly undistinguished CTTS prior 
to its 2010 eruption. 

3. 2009 AND 2010 OBSERVATIONS 
3.1. Optical Photometry and an Independent Discovery 

During the 2009 and 2010 observing seasons, PTF ob- 
tained red optical images of the North America/Pelican 
Nebula star-forming region with a typical 5-day cadence. 
These observations were conducted with the main PTF 
Survey Camera, the former CFHT12K mosaic camera 
now extensively re-engineered and mounted on the 48-in 
Samuel Oschin telescope at Palomar Observatory (here- 
after P48). The camera is a mosaic of 12 CCDs (one of 
which is not functional), covering a 7.8 square degree field 
of view with 1" sampling; typical conditions at Palomar 
Observatory produce 2.0" full width at h alf-maximum 
intensity (FWHM) images ( |Law et al.||2009[ ). The i? PT F 
filter, a Mould R band, is similar to the Sloan Digital Sky 
Survey (SDSS) r band in shape but is shifted redward 
by ~ 27 A and is ~ 20 A wider. The typical 5a limiting 
magnitudes are ttlr ~ 20.6 (AB) in 60 s exposures. Rep- 
resentative images of LkHa 188-G4/PTF lOqpf from the 
2009 and 2010 observing seasons are shown in Figure [2j 

Transient sources are detected in the PTF monitoring 
data by means of automated reduction pipelines, includ- 
ing a nearly real-time image-subtraction pipeline hosted 
at Lawrence Berkeley National Laboratory (LBNL). 
Well-detected sources in the difference images are scored 
(using a human-trained, machine-based classifier) for 
their likelihood of being truly astrophysical in nature or 
of spurious origin. Variable sources with larger likeli- 
hoods of being nonspurious are passed to an automatic 
source classifier at UC Berkeley ( Oarical), which com- 
bines PTF measurements with all other available infor- 
mation (e.g., SIMBAD identifications, 2MASS photom- 
etry, etc.) to provide probabilistic classifications of PTF 



4 



Miller et. al. 




• • • 



50 100 150 200 250 300 350 400 450 500 
MJD - 55000 

Fig. 3. — P48 light curve showing the outburst of LkHa 188- 
G4/PTF lOqpf. During the 2009 season PTF lOqpf was variable 
with typical Am PS 0.2 mag. The source dramatically increased the 
rate at which it was brightening around 2010 July 13 (day ~391 
in the figure). P48 saturates on sources brighter than mn Fi „ Ri 
13.5 mag. The shaded region shows all data where PTF lOqpf is 
saturated; these data represent lower limits to the true brightness 
of PTF lOqpf during those epochs. 



detections (Bloom et al., 2011, in prep.). These initial 
classifications are made available to PTF collaboration 
members via the PTF Follow-up Marshal, which enables 
visual inspection of current and reference images, precur- 
sor PTF light curves, and any subsequent spectroscopy. 

The source reported here was detected by the 
PTF pipeline and automatically assigned the name 
PTF lOqpf, following internal PTF naming conven- 
tions. Differential photometric measurements are made 
via aperture photometry, with the aperture size scaled 
to the seeing in each image, using SExtractor (Bertin 
& Arnouts||1996|), and absolute calibration is relative to 
BT 



the USJNO-B1 catalog ( |Monet et al. | 2003 [) with an uncer 

ine 



tainty of ~ 0.15 mag. CCD 3 on the P48 camera satu- 
rates at 55,000 counts, roughly corresponding to sources 
with R < 13.5 mag, given the standard 60 s exposure 
time. Thus, in all of our observations following those on 
2010 Aug. 14, LkHa 188-G4 is saturated and our mea- 
surements represent a lower limit to the true brightness 
of the source. These limits are reflected in the full P48 
light curve of LkHa 188-G4, shown in Figure [3j and the 
photometry is reported in Table[5| We do not include the 
systematic ~ 0.15 mag calibration uncertainty in this ta- 
ble because this uncertainty represents a constant offset 
that would be applied to all the data in the same way. 
The PTF also makes use of the robotic Palomar 60- 

to obtain multi- 



2006) 



in telescope (P60; Cenko et al. 

filter photometry for source verification and classifica- 
tion purposes. PTF lOqpf was observed with the P60 on 
2010 Oct. 14; Figure [4] presents a three-color image con- 
structed from the P60 riz frames. Table [2] contains the 
individual photometric measurements. The absolute ze- 
rop oint calibration was done relativ e to SDSS photome- 
try ( Adelman-McCarthy et al. 2008 ) based on other fields 
observed by the P60 during the same night with the same 
filter. The uncertainty in the zeropoint is computed as 
the standard deviation of the zeropoints in all other SDSS 
fields. 



Fig. 4. — False-color riz image of PTF lOqpf during the 2010 
eruption. The image is 2.2' on a side with north up and east to the 
left. PTF lOqpf is highlighted with crosshairs. Notice the bright 
reflection nebula to the southwest of PTF lOqpf; this nebula was 
not present prior to the 2010 eruption. 



3.2. Near-Infrared Photometry 

Near-infrared observations of PTF lOqpf were con- 
ducted with the 1.3 m Peters Automated Infrared Im ag- 
ing Telescope (PAIRITEL; 



Bloom et al 



Mt. Hopkins, AZ, starting on 2010 Sep. 27. PAIRITEL is 
a roboticized system using the former 2MASS southern 
hemisphere survey camera that employs two dichroics to 
observe simultaneously in the J, H, and K s bands. Ob- 
servations were scheduled and executed via a robotic sys- 
tem. PAIRITEL is operated in a fixed observing mode 
in which 7.8 s double-correlated images are created from 
the difference of a 7.851 s and a 51 ms in tegration taken 
in rapid succession (see Blake et al.|20 08 ) . The standard 
observing procedure involves taking three image pairs 
prior to dithering the telescope. 

The raw data from these images are reduced us- 
ing standard IR reduction methods via PAIRITEL 
PIPELINE III and the flux for all sources is measured 



via aperture photometry using SExtractor (Bertin & 
Arnouts||1996[) , calibrated against 2MASS. In the H and 
K s bands, PTF lOqpf saturates the 7.851 s frames; how- 
ever, PIPELINE III produces "short-frame" mosa ics con- 
sisting of re duced, stacked 51 ms images (see also Bloom 



et al. 20091). The "short-frame" mosaics contain >10 
bright 2MASS stars which we use to properly calibrate 
photometric measurements of PTF lOqpf in these im- 
ages. PAIRITEL has a systematic uncertainty of ^0.02- 



0.0 3 mag in each of the J, H , and K s bands (see Blake et 



al. 2008 



Perley et al.||2010 ) , which, in the case oFLkHa 



188-C4, is larger than the statistical error in all three 
bands. Thus, we add a systematic error of 0.03 mag in 
quadrature with the statistical uncertainty to determine 
the total uncertainty in each band. PAIRITEL photom- 
etry is reported in Table [3j and a J, H, and K s band 
false-color image of PTF lOqpf is shown in Figure [5] 

3.3. Optical Spectroscopy 

3.3.1. Lick Spectroscopy 

Low-resolution spectra of PTF lOqpf were obtained on 
2010 Sep. 16 and 2010 Nov. 2 with the Kast spectro- 
graph on the Lick 3-m Shane telescope (Miller & Stone 
1993). The spectra were r educed and calibrated using 
standard procedures (e.g., |Matheson et al. 2000). For 
the Sep. 16 spectrum, flux calibration for the red arm 
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flat-field, wav elength calibrated, and optimally extracted 
(Home 11986). The spectra were flux-calibrated and cor- 
rected for telluric absorption using observations of an 
AOV star at a similar airmass with the IDL-based code 



Fig. 5.— False-color JHK a image of PTF lOqpf during the 2010 
eruption. The NIR reflection nebula is less prominent than the 
optical nebula, but still visible to the southwest of PTF lOqpf. The 
image is 5' on a side with north up and east to the left. PTF lOqpf 
is highlighted with crosshairs. 

of the Kast spectrograph was performed using a stan- 
dard star observed at high airmass, whereas PTF lOqpf 
was observed at low airmass; thus, the absolute flux cal- 
ibration is somewhat uncertain in the red portion of the 
optical spectrum. Furthermore, the noisy features near 
-9300 A arc likely the result of an imperfect telluric cor- 
rection and not astrophysical. Since conditions were not 
photometric on the night of Nov. 2, the absolute flux cal- 
ibration is again somewhat uncertain. Nevertheless, on 
both nights we observed PTF lOqpf with the slit placed 
at the parallactic angle, so the relative spectral shapes 
should be accurate. The full low-resolution spectra of 
PTF lOqpf are shown in Figure [6j while all of our spec- 
troscopic observations are logged in Table |4j 

3.3.2. High- Resolution Spectroscopy 

PTF lOqpf was observed on 2010 Sep. 25 with the 
Keck I telescope and HIRES spectrometer (Vogt et al. 
1994 1 using the red cross-disperser a nd the standard set 



tings of the Californ ia Planet Survey ( |Howard et aLpOlO 
Johnson et al. 12010). The C2 decker was employed, pro 



jecting to 078b on the sky and providing a resolution 
R = A/AA = 55,000 at 5500 A. We took advantage 
of the standard CPS reduction pipeline, which includes 
fiat-fielding, scattered-light subtraction, order tracing, 
cosmic ray rejection, and spectrum extraction. After 
extraction, each order is sky-subtracted and summed 
in the cross-dispersion direction to form the final one- 
dimensional spectrum. The exposure time of 560 s re- 
sulted in a signal-to-noise ratio (S/N) of - 20 at 7400 A. 

3.3.3. Near- Infrared Spectroscopy 

A NIR spectrum of PTF lOqpf was obtained with the 
Triple Spec spectrograph o n the Palomar 5-m Hale tele- 
scope QHerter et aL]|2008[ ) on 2010 Sep. 23. TripleSpec 
has no moving parts and simultaneously acquires 5 cross- 
dispersed orders covering 1.0-2.4 (im at R ~ 2700. 

The spectra were reduced with an IDL-based data re- 
duction pipeline developed by P. Muirhead 12 . To facil- 
itate the subtraction of the sky and background emis- 
sion of the total signal, the observations were obtained 
using an ABBA dither pattern along the slit. Each sky- 
subtracted exposure was then divided by a normalized 

12 http://www.astro.cornell.cdu/~muirhcad/#Downloads. 



xtellcor .general by Vacca, Cushing, & Rayner (2003). 
The full TripleSpec spectrum is shown in Figure [7j 

4. ANALYSIS 

4.1. Image Morphology and Photometric Analysis 

The large-amplitude eruption (—4 mag) of LkHa 188- 
G4/PTF lOqpf qualitatively matches the defining char- 
acteristic of FU Orionis eruptions. Such an eruption, 
however, is not sufficient to prove that a YSO is in an 
FU Orionis-like state. PTF lOnvg, for example, bright- 
ened by —6 mag over the course of a year, yet shows lit- 
tle resemblance to FU Ori-like stars (Covey et al. 2011). 
Spectroscopy over a broad wavelength range (see below) 
is required to conclusively classify a source as an FU 
Orionis-like star. 

Over the course of our P48 observations, PTF lOqpf 
brightened by ~1 mag from early August 2009 to early 
August 2010, before suddenly brightening by ~3 mag 
over a ~2 month period (see Figure p|. Between the start 
of our observations on 2009 Aug. lJand 2010 Jul. 9, the 
initial slow rising phase, PTF lOqpf rose by an average 
of -0.003 mag d" 1 , while between 2010 Jul. 9 and 2010 
Aug. 11, during the rapid rising phase, PTF lOqpf rose 
by an average of —0.07 mag d _1 . 

After 2010 Aug. 11, PTF lOqpf is saturated in the P48 
images and we cann ot constrain the rat e at wh ich it is 
increasing; however, ISemkov & Peneva (2010b) report 
that it reached peak brightness on 2010 Aug. 24 with 
R = 12.79 mag. Despite the rather rapid rise at the be- 
ginning of the PTF lOqpf optical outburst, the overall 
light curve is remarkably smooth. The P48 light curve 
for PTF lOqpf shows fluctuations at the ~0.2 mag level 
about both the pre-outburst slow rise and the sharp rise 
of the outburst. The long-term trend with low scatter 
observed prior to 2010 July is somewhat unusual for a 
standard CTTS, whose year-to-year variability is typi- 
cally significantly sma ller than the intra-year variability 
( |Gr'ankin et al |[2007[ ). Indeed, the IPHAS photometry 
of LkHa 188-0-4 captures variability that is more char- 
acteristic of CTTSs: variability at the ~0.3 mag level 
was observed over the course of 1 month in 2003, with 
a follow-up observation two years later falling within the 
same range of magnitudes as observed in 2003. This 
suggests that the slow rise LkHa 188-G4 demonstrated 
in 2009 was a prelude to the 2010 outburst, rather than 
simply standard CTTS variability. 

A new reflection nebula, which was not present in 2009, 
can now clearly be seen around PTF lOqpf (Figure [2J. 
Our regular observations over the course of the rise of 
PTF lOqpf allow us to constrain the appearance of the 
nebula, si milar to the analysis of V1647 Ori/McNeil's 



Nebula by Bricefio et al. ( 2004 ) . In the P48 subtraction 
images produced as part of the PTF reduction pipeline, 
the residual is asymmetric about the point-source cen- 
troid on 2010 Aug. 2. The dust responsible for the re- 
flection nebula is offset to the southwest of PTF lOqpf, 
and the asymmetry in the subtraction image is the re- 
sult of this dust scattering light into our line of sight. 
We measure R « 14 mag for PTF lOqpf on 2010 Aug. 
2. The nebula may be present as early as 2010 July 19, 
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Wavelength (A) 

Fig. 6. — Low-resolution spectra of PTF lOqpf. Prominent absorption lines are marked. Ho is the only line seen in emission, and it 
appears to have decreased in strength between our two spectra. Residuals left following an imperfect removal of telluric absorption have 
been marked with a ©, while the noisy features around ~9300 A are also likely telluric (see text). The spectrum from 2010 Nov. 2 has 
been shifted down by 0.5 for clarity. 



however: the faint nebula is difficult to distinguish from 
the wings of the central star's point-spread function. The 
optical reflection nebula (Figure [4]) has an apparent size 
of ~40", which at a distance of ~600 pc corresponds to 
~24,000 AU. 

In the NIR we have less information than in the optical. 
Our observations suggest a flat light curve in the J, H, 
and K s bands starting on 2010 Sep. 27. Our NIR photo- 
me tric measurements are consistent with those presented 
by Leoni et al. (2010), and we confirm that PTF lOqpf 
has brightened by 3.25, 3.07, and 2.81 mag in the J, H, 
and K s bands, respectively, since 2MASS observations 
were obtaine d on 2000 Jun. 10. T he NIR measurements 
presented by Cohen & Kuhi ( 1979 ) are ~1.1 mag brighter 
than those of 2M ASS, wh ich may ~b e the result of a larger 
beam size for the Cohen & Kuhi (19791 measurements. 
FU Orionis star s are frequently a ssociat ed with NIR re- 



flection nebulae ( Connelley, Reipurth, fc Tokunaga|2007 



Connelley fc G reene 2010}, and PTF lOqpf is no excep 
tion (see Figure |5|. The NIR nebula is similar in size to 
the nebula seen in the optical. 

The outburst SED is shown in Figure [l] and it clearly 
demonstrates that PTF lOqpf is bluer in the optical/NIR 
than it was prior to the 2010 outburst. We have limited 
spectral coverage during the outburst; in particular, we 
lack any information on how the emission has changed in 
the MIR. Direct integration of the optical/NIR outburst 
SED results in a lower limit to the luminosity of at least 
4 Lq. Assuming the same bolometric correction to the 
H band before and after the outburst results in an out- 
burst luminosity of ^outburst ~ 12 Lq. This luminosity 
places PTF lOqpf at the bottom of the luminosity scale 
for FU Ori outbursts, and more t han an order of mag- 



4.2. Spectroscopic Analysis 

Optical spectra of FU Orionis variables share a few 
defining properties: a late F or G supergiant photo- 
sphere, an absorption-dominated spectrum with only Ha 
typically seen in emission, and broad P Cygni absorption 
with representative velocities of a few hundred km s . 
Figure [6] shows that PTF lOqpf exhibits these general 
properties. Figure [9] illustrates our HIRES spectra cov- 
ering three Balmer transitions. Each of these lines ex- 
hibits a strong P Cygni absorption profile, with absorp- 
tion minima around —40 km s _1 . Ha is the only line 
having an emission component, and this emission moves 
the absorption minimum slightly blueward of the other 
Balmer lines. The absorption extends to < —200 km s _1 , 
suggesting in the classical interpretation that PTF lOqpf 
is driving a strong outflow (see also Figure [9]). 

In the discussion below, all velocities havebeen trans- 
formed to the local standard of rest (LSR) . To our knowl- 
edge, a precise radial velocity measurement for LkHa 
188-G4 has not been made. Therefore, we adopt a rest 
velocity of VLsr ~ 1-6 km s -1 , consistent with 13 CO ob- 



servations of the surrounding molecular gas ( Dobashi et 
al p994| ). 



nitud e fainter than FU Ori itself (Hartmann & Kenyon 
19961. 



4.2.1. Low-Resolution Optical Spectroscopy 

At low resolution, PTF lOqpf exhibits little evolution 
between our two Kast spectra taken 40 d apart (see Fig- 
ure [6]). The strength of the absorption features remains 
constant to within the uncertainties, while the only ma- 
jor change is in the emission flux from Ha. During this 
40 d period the Ha flux has decreased from ~ 2.8 x 10 30 
erg s" 1 to ~ 1.5 x 10 30 erg s _1 . Observing at the paral- 
lactic angle allows us to accurately measure the change 
in the relative color during the 40 d gap between ob- 
servations. Given the typical uncertainties associated 
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Fig. 7. — TripleSpec NIR spectrum of PTF lOqpf showing from top to bottom the YJ, H, and i\-band specta. Prominent absorption 
features in the outburst spectrum are identified, w hile the spectrum shows no clear em ission lines. For comparison the spectra of M5 lb II 
and K5 lb supergiants, and an M0 V dwarf (from [Rayner, Cushing, fc V acca 2009), are also shown. The muted features of PTF lOqpf 
relative to the M0 V star shows that the NIR emission is trom a low surface gravity environment. PTF lOqpf exhibits strong 2.3 fim CO 
absorption and ~1.4 and 1.9 fim H2O absorption, similar to other FU Orionis-like sources. 



with spectre-photometric measurements with the Kast 
spectrograph (Silverman et al., 2011, in prep.), we find 
that PTF lOqpf does not show significant evidence for a 
change in color between our two observations. 

To constrain the optical spectral type of PTF lOqpf, 
we use a least-squares minimization procedure to find 
the best match between the low-resolution Kast spec- 
trum and the low-resolution stellar spect ra provided in 
the Pickles Stellar Spectral Flux Library ( |Pickles|[l998] ). 
The spectral type of FU Orionis v ariables is kn own to 
vary as a function of wavelength (Herbig 1977); thus, 
we divide the spectra into two sections, 4000-6000 A 
and 6000-9000 A. Below 4000 A, the S/N is low in our 
Kast spectrum, while above ^9000 A there is an uncer- 
tain telluric correction which may add artificial features. 
Next we redden each of the Pickles star spectra in in- 
crements of 0.05 mag from E(B-V) — 0.00 to 2.25 mag 
with Ry =3.1. We mask the region around Ha, since 
emission lines are not expected in normal stellar spectra. 
We then perform a least-squares fit of the reddened li- 
brary spectra with our low-resolution spectrum of LkHa 



188-G4, where the only free parameter is an overall nor- 
malization constant. 

We illustrate the results of this procedure in Fig- 
ure |5J where we include the low-resolution spectrum 
of PTF lOqpf along with the four best matches from 
our least-squares minimization. The stellar spectra are 
shown in the order of the quality of the fit, with the best 
fit on top and the fourth-best fit on the bottom. The 
adopted reddening has little physical meaning given that 
the observed flux is not emanating from a stellar photo- 
sphere having the standard single-temperature SED; this 
portion of the procedure merely serves as an automated 
process to vary the slope of the continuum, which is nec- 
essary given that PTF lOqpf has a heavily extinguished 
spectrum. 

The Pickles library is incomplete — it does not in- 
clude every luminosity class for every spectral type. For 
this reason we emphasize that our procedure is only il- 
lustrative; it does not provide a definitive spectral type 
for PTF lOqpf in cither the blue or red portion of the 
optical spectrum. Nevertheless, a few trends are readily 
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Fig. 8. — Low-resolution optical spectra of PTF lOqpf compared to stellar spectra from the Pickles library (Pickles 1998). The comparison 
spectra are selected as the best matches to PTF lOqpf following a least-squares minimization procedure (see text). The comparison spectra 
are shown in order of the goodness of fit, with the best match on top. The reddening values were used as a proxy for removing the overall 
continuum shape. 



apparent from Figure [8] In the blue portion of the op- 
tical, PTF lOqpf most closely resembles the spectra of 
G supergiants and G giants. In the red portion of the 
optical, the spectrum resembles that of G and K stars, 
with no real hint regarding the luminosity class. The 
adopted values of the reddening are significantly different 
between the blue and the red, but again we remind the 
reader that this reddening is only a device with which 
to match the intrinsic shape of the continuum. Typi- 
cally, the Mg I A5172 triplet is used to determine the 
luminosity class of G-type stars, but our high-resolution 
spectra show that these lines are tracing an outflow (see 
Figure [9]) and thus they are unsuitable for providing a lu- 
minosity class. Despite this, the procedure does demon- 
strate that PTF lOqpf is similar to G giants/supergiants 
in the blue portion of the optical, consis tent with ob- 
servations of other FU Orionis variables (Hartmann & 
Kenyon||l996| . Additionally, as we will show below, the 
spectrum shifts to later types (K stars in this case) with 
increasing wavelength, as is also observed in FU Orionis 
variables (Herbig 1977). 

4.2.2. High- Resolution Optical Spectra 

The HIRES spectrum of PTF lOqpf is dominated by 
broadened stellar-like absorption features. The spectral 
lines are consistent with those of a comparison G2 I 
star and inconsistent with spectral types later than early 
K. We note, however, that the diagnostics of G vs. K 
types are subtle and would become blurred with the large 
broadening of the PTF lOqpf lines. The line broadening 
is similar to, but perhaps slightly less than, the 60-100 
km s _1 broadening of V1057 Cyg, and significantly more 
than the 30-40 km s _1 of V1515 Cyg in which some hints 
of a K-type spectrum can be seen at redder optical wave- 
lengths. A detailed comparison of PTF lOqpf to these 
two canonical FU Ori objects reveals that the spectrum 
of PTF lOqpf is similar to both in terms of the relative 
strength of the common absorption features. However, 
we note that each object has unique and notable spec- 
troscopic attributes. 

Several lines have clear kinematic signatures of winds, 
specifically Ca II H & K, Fe II A5018, the Mg I AA5167, 



5172, 5183 triplet, the Na I D lines, Ha, Hy, H<5 (H/3 falls 
in a gap between the orders in our HIRES setting), and 
the K I A7699 line (its doublet partner at 7665 A also 
falls in a gap between orders). The morphology of these 
lines is illustrated in Figure |9j Notably, the blueshifted 
absorption troughs are narrower in our early-stage spec- 
trum of PTF lOqpf than exhibited in modern spectra of 
V1515 Cyg and V1057 Cyg, by at least a factor of two. 
Altho ugh LkHa 188-G4 ex hibited [O I] A6300 emission 
in the|Cohen & Kuhi (1979) spectrum, it is not apparent 



post-outburst, even at high spectral resolution, nor do 
we see other forbidden species in emission (such as [N II] 
and [S II]) typically associated with outflows in YSOs. 

The other line worth noting in PTF lOqpf for its kine- 
matic signatures is Li I A6707. Figure [TO] shows the spec- 
tral region containing this line, in comparison to V1515 
Cyg and V1057 Cyg as well the G2 I standard. The Li 
I absorption profile has multiple components including a 
narrow slightly blueshifted component arising from low- 
velocity gas and a broader absorption extending to ±60 
km s _1 that is reminiscent of the flat-bottomed profiles 
seen in other metallic lines in the same spectral region. A 
similar absorption profile is also seen in Ba II (Figure [9]). 

Figure [Tl] illustrates the flat-bottomed nature of many 
absorption lines in PTF lO qpf. The profiles are re minis- 
cent of those illustrated by |Petrov fe He rbig ( 2008]) for 
FU Ori. When the flat-bottomed lines in PTF lOqpf are 
cross-correlated against a normal G2 I stellar template, 
the resulting cross-correlation function is broad (~120 
km s" 1 ) with several different peaks. In some spectral re- 
gions a central peak can be fit by a Gaussian of width 30 
km s^ 1 , and there are two secondary peaks near the full 
width of the correlation function, each having Gaussian 
width ~20 km s _1 . In other spectral regions, both blue 
and red, the correlation function has no central peak, but 
only the double peaks near the limb. To measure the 
goodness of the cross correlation we adopt the Tonry- 
Davis R value (Tonry & Davis 19791, which is a ratio 
of the power in the cross-correlation peak to a measure 
of the antisymmetry in the correlation. Large R values 
indicate good cross correlations. Using a G2 supergiant 
as a template, the Tonry-Davis R value for PTF lOqpf 



A New FU Orionis Outburst From a CTTS 



9 




Fig. 9. — Keck HIRES high-resolution profiles of several prominent features showing strong blueshifted absorption, indicative of outflowing 
material. Top left: The Balmer series for PTF lOqpf. H/3 fell in a gap between echelle orders for our observing setup, and thus is not 
displayed. PTF lOqpf exhibits strong P Cygni absorption in the other Balmer lines with outflow speeds of ~100-300 km s _1 . Top center. 
Profiles of the Mg I triplet at 5167, 5172, and 5183 A. The lines show strong asymmetries, tracing the kinematics of the wind from 
PTF lOqpf. Fe II A5169 is present at about +60 km s relative to Mg I A5167. Top right: Absorption profiles of the Na I D doublet 
showing similar kinematics as the Mg I lines. Bottom left: Absorption profiles of Ca II H &; K. He is also present, offset from Ca H by 
~100 km s _1 . Bottom center: Absorption profiles of Fe II A5018 and K I A6989. Bottom right: Absorption profiles of Ba II AA4554, 5853, 
and 6141. These lines show narrow absorption, FWHM ps 20 km s _1 , superposed on broader, shallower profiles extending ± ~6 km s 
from the line c enter. For clarity, all spectra have been smoothed with a Savitzky-Golay filter, which is similar to a running mean ( |Savitzkyj 
|fc Golay|1964} . Velocities are given relative to the LSR, and the vertical grey line indicates the rest velocity of PTF lOqpf, Vlsr ~ 1-6 



km s ^see text). 
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Fig. 10.— Keck HIRES high-resolution spectra of the Li I A6706 
and Ca I A6717 region of PTF lOqpf. For comparison we also show 
high-resolution spectra of V1057 Cyg, V1515 Cyg, and a G2 lb 
spectral standard, all from our unpublished spectral database. In 
particular, notice the strong Li absorption present in PTF lOqpf 
and the two canonical FU Orionis objects. The Li line in both 
V1057 Cyg and PTF lOqpf traces an unusual kinematic structure, 
possibly similar to that seen in Ba II (see Figure [9]) 



is 4-6 (compared to R > 5 to 15 for V1057 Cyg, which 
shows broadening and flat-bottomed structure similar to 
PTF lOqpf in some lines, and R > 20 for V1515 Cyg, 
which does not). The cross-correlation results are sim- 
ilar when using V1515 Cyg as a template, rather than 
the G2 I standard. Employing V1057 Cyg as a template 
results in a broad and centrally peaked cross-correlation 
function having width 120-130 km s _1 . 

4.2.3. Near-Infrared Spectroscopy 

NIR spectra of FU Orionis variables resemble the spec- 
tra of K and M giants/supergiant s, in contrast to their 



earlie r spectral type in the optical ( Hartmann & Kenyon 
1996|. The most prominent features in the JN1R spectra 



of F IT Orionis stars are the strong absorption bands from 
CO at 2.3 /iin and the H 2 vibration-rotation bands at 
1.38 /im and 1.87 /im. These absorption bands are clearly 
present in the outburst spectrum of PTF lOqpf; see Fig- 
ure [71 which compares spectra from |Rayner, Cushing, fc| 



Vacca 



2009| ) of two supergiants (K5 lb and M5 Ib-II) 



and a dwarf (M0 V) with PTF lOqpf. There are many 
prominent absorption features seen in the M0 V spec- 
trum, which are muted or absent in PTF lOqpf and the 
two supergiant spectra. This is indicative of a low surface 
gravity environment for PTF lOqpf. 
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Fig. 11. — Several Fe I line profiles illustrating the flat-bottomed 
nature typical of the other metallic lines. The Fe I transitions at 
AA6393, 6593, and 6678 have been shifted to the same velocity, 
while Fe I A6400 is offset by ~300 km s" 1 relative to A6393. These 
lines all exhibit a prof ile with a "flat" bott om similar to those 
observed in FU Orionis (IPetrov & Herbig|2008|l. For clarity, spectra 
have been smoothed with a Savitzky-Golay niter, which is similar 
to a running mean. The vertical grey line is the same as in Figure|9] 



We have attempted to derive a crude NIR spectral type 
for PTF lOqpl. In addition to its anomalously strong CO 
features, PTF lOqpf also possesses prominent H 2 ab- 
sorption in the H and K bands; the relative strengths 
of the H and K-hanA absorptions are unusual, with 
PTF lOqpf separated from both the d warf and giant loc i 
for those features (see, e.g., Figure 4 of Covey et al. 2010 ). 
We therefore attempted to characterize the JN1R spectral 
type of PTF lOqpf based on the strengths of several H- 
band absorption features: Mg I (1.51 (an), K I (1.518 
/mi), Al I (1.67 /mi), and a blend at 1.70 /mi. We com- 
pared the strengths of these features to th ose measured 
for a s uite of dwarf and giant templates by |Covey et al. 
(2010). While this method is sufficient to derive only a 
crude NIR spectral type for PTF lOqpf, we find that the 
source's outburst spectrum most closely resembles that 
of K7 III — M2 III giants (see Figure n\ for some com- 
parisons). It is important to note that the classification 



scheme from Covey et al. ( |2010 ) does not include su- 
pergiants; thus, we conclude that PTF lOqpf has a low 
surface gravity, but we cannot discriminate among lumi- 
nosity classes I— III. 

While the strength of molecular water absorption is 
greater than that normally observed in giants, this find- 
ing is con sistent with observ ations of other FU Orio- 
nis stars. |Sato et al.| (|1992|) argue that the unusually 



strong H 2 absorption seen in FU Orionis stars can be 
explained if water that has condensed onto dust grains in 
the midplane of the disk becomes vaporized during the 
heating of the disk associated with an FU Orionis erup- 
tion. Once this evaporated water is flung vertically out 
of the disk midplane, it will imprint a strong absorption 
si gnature in the NIR spectrum . 

Connelley & Greene (2010) show that the equivalent 
width (EW) of CO absorption at 2.3 /mi can be used to 
identify FU Orionis-like stars when compared to the EW 
of Na (2.206 and 2.209 /mi) and Ca (2.263 and 2.266 
/mi). To compare PTF l Oqpf to the stars studied by 
Connelley & Greene (2010), we measure the EW of Na, 



Ca, and CO fo llowing the prescription of|Chalabaev fe 



Maillard fl!983|) with the correction provided by |Voll- 



maim & Eversberg|(|2006|). We measure Na and Ca EWs 



of 1.56 ± 0.05 and 1.07 ± 0.06 A, respectively, while the 
EW of CO is 35.6 ± 0.3 A. These uncertainties reflect 
the statistical errors in our measurements, and likely un- 
derestimate the systematic effects associated with these 
measurements (e.g., the precise placement of the contin- 
uum). Nevertheless, the systematic uncertainties are not 
large enough to alter the general conclusions drawn be- 
low. The CO absorption is well in excess of what would 
be expected for both late-type dwarfs and giants based 
on the strength of the Na and Ca absorption; however, 
this excess CO absorption is consistent with other FU 
Orionis-like stars in the sample of Connelley & Greene 
(see their Figure 4). 

Prominent H and He I A10830 absorption lines are also 
visible in the NIR spectra (see Figure [7]). However, H 
and He absorption cannot originate in a simple, single- 
temperature, cool (late K-M) atmosphere. Thus, these 
lines, though unresolved in our spectra, likely trace an 
outflow, similar to the Balmer series in the optical. He I 
A10830 may be a hot-wind indicato r, but high-resolution 



spectra are needed to confirm th is ( Edwards et al.|[2006 
Kwan, Edwards fe Fischer]|2007[ ) . 



5. DISCUSSION AND CONCLUSIONS 

We have shown that the 2010 eruption of LkHa 188-G4 
places it within the small family of FU Orionis-like stars 
having observed outbursts. LkHa 188-G4/PTF lOqpf 
exhibits the defining observable characteristics of the 
FU O rionis class as described by |Hartmann fe Kenyon| 
(fl996|: 



• it exhibited a > 4 mag increase in its optical light 
curve; 

• it is associated with a star-forming region and 
shows a bright reflection nebula following its erup- 
tion; 

• optical spectra of the outburst are consistent with 
a G-type giant /supergiant; 

• during outburst, the Balmer lines exhibit strong P 
Cygni profiles, with only Ha showing an emission 
component, and absorption extending to about 
-200 km s" 1 ; and 

• the NIR spectra during the eruption resemble those 
of M-type giants/supergiants, with strong absorp- 
tion seen in the H2O and CO bands. 

In sharp contrast to the class of YSO outbursts referred 
to as EX Lupi-like or V1647 Ori-like, PTF lOqpf docs 
not exhibit strong optical/NIR emission lines during out- 
burst. 

In the IR, unlike typical M and K-type supergiants, 
PTF lOqpf shows prominent absorption lines of H I (Br7, 
Pa/3, Pa7, Pa,8) and He I 1.083 /tm. In the optical the H 
I lines trace an outflow, with clear P Cygni absorption 
profiles observed in our HIRES spectrum. Our spectral 
resolution is insufficient to determine whether a similar 
P Cygni profile is observed in the NIR; however, H and 
He absorption cannot originate in an M-giant/supergiant 
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photosphere. Thus, the H and He lines are likely the 
result of inner-disk region accretion/wind physics. 

LkHa 188-G4 stands out among the currently known 
group of FU Orionis-like stars in that there are pre- 
outburst observations with detections from the optical 
to the MIR, making LkHa 188-G4 the first FU Orionis- 
like star with a well-sampled SED prior to eruption. This 
SED shows that LkHa 188-G4 was a Cla ss II YSO, w hich 
are commonly associated with CTTSs (Lada 1987). A 



pre-outburst optical spectrum reveals a late-type star, 
K7-M0, with prominent Ha and H/3 seen in emission, 
again consistent with a CTTS. Multiple epochs of optical 
imaging prior to the outburst reveal that LkHa 188-G4 
exhibited stochastic variability at a level consistent with 
a CTTS. 

One of the most important elements in the study 
of PTF lOqpf/LkHa 188-G4 is that it is associated 
with a known, optically revealed and previously stud- 
ied YSO. The clear identification of LkHa 188-G4 as a 
CTTS prior to eruption solidifies the interpretation of 
FU Ori events as enhanced accretion and outflow likely 
associated wit h disk-accretion instabilities ( Hartmann & 
Kenyon 19961. This outburst's association with a pre- 
viously classified CTTS source shows that FU Ori erup- 
tions are not strictly limited to the most highly embedded 
Class 0/1 phases of pre-main sequence evolution. Indeed, 
the strong winds associated with these events may play 
an important role in disrupting a YSO's circumstellar 
envelope. Such a disruption is thought to be a key com- 
ponent of the Class I /II transition. 

Continued photometric and spectroscopic monitoring 
of PTF lOqpf/LkHa 188-G4 will fully elucidate its na- 
ture. As more and more synoptic surveys come on- 
line, the rich legacy of IR observations made by sev- 
eral ground- and space-based observatories (e.g., 2MASS, 
IRAS, MSX, Spitzer, AKARI, WISE) over the past two 
decades will enable the discovery of additional FU Ono- 
nis outbursts with well characterized pre-outburst SEDs. 
These discoveries will in turn allow us to determine if the 
case of LkHa 188-G4 is an outlier or the norm. A bet- 
ter understanding of whether FU Orionis-like eruptions 
typically occur during the Class I or Class II phase of pre- 
main sequence evolution will allow us to place stronger 
constraints on the triggering mechanism and occurrence 
rates for these events. 
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TABLE 1 

IPHAS Observations of LkHo 188-G4. 



UT date 



2003/10/15 
2003/10/18 
2003/11/02 
2003/11/14 
2005/11/19 



r mag 
(AB) 
17.300±0.010 
17.348±0.011 
17.283±0.017 
17.074±0.009 
17.075±0.007 



i mag 

(AB) 
15.677±0.006 
15.744±0.006 
15.608±0.009 
15.570±0.005 
15.671±0.006 



Ha mag 
(AB) 
16.023±0.007 
16.063±0.006 
15.975±0.009 
15.724±0.006 
15.768±0.008 



TABLE 2 
P60 Observations of PTF IOqpf. 



date 


filter 


mag 


°"mag 


(MJD) 




(AB) 




55483.127 


r 


12.948 


0.039 


55483.127 


i 


12.207 


0.027 


55483.128 


z 


11.623 


0.041 



TABLE 3 

PAIRITEL Observations of PTF IOqpf. 

*mid a J mag H mag K s mag 

(MJD) (Vega) (Vega) (Vega) 

55466.137 10.04 ± 0.03 9.14 ± 0.03 8.65 ± 0.03 

55468.145 9.99 ± 0.03 9.06 ± 0.04 8.64 ± 0.04 

55469.148 10.04 ± 0.03 9.07 ± 0.03 8.70 ± 0.03 

55479.109 10.11 ± 0.03 9.11 ± 0.03 8.63 ± 0.04 

55504.164 10.20 ± 0.03 9.24 ± 0.03 8.74 ± 0.03 

55513.195 10.29 ± 0.03 9.34 ± 0.03 8.79 ± 0.03 

55518.168 10.32 ± 0.03 9.34 ± 0.04 8.84 ± 0.03 

55527.117 10.35 ± 0.03 9.40 ± 0.03 8.89 ± 0.03 

55531.145 10.40 ± 0.03 9.44 ± 0.03 8.97 ± 0.03 

55543.066 10.45 ± 0.03 9.48 ± 0.04 9.06 ± 0.04 

"Midpoint between the first and last exposures in a single 
stacked image. 



TABLE 4 

Log of spectroscopic observations. 



UT Date 


Telescope/ 


Exposure 


A 




Instrument 


(s) 


(A) 


2010-09-16.24 


Shane 3-m/Kast 


600 


3400-10300 


2010-09-23.24 


Hale 5-m/TripleSpec 


720 


9400-24,600 


2010-09-25.30 


Keck I 10-m/HIRES 


560 


3640-7990 


2010-11-02.21 


Shane 3-m/Kast 


1200 


3400-10300 
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TABLE 5 

P48 i?-BAND Photometry of PTF IOqpf 



date mag o"mag 
(MJD) 



5 5 5 6 


221 


17 


023 





017 


5 5 5 6 


295 


17 


098 





019 


55061 


383 


16 


843 





012 


55062 


327 


16 


880 





020 


55064 


286 


16 


950 





018 


55064 


331 


16 


818 





021 


55067 


393 


16 


803 





031 


55080 


361 


17 


079 





036 


55093 


318 


16 


933 





016 


55094 


252 


16 


777 





018 


55107 


222 


16 


981 





026 


55107 


266 


16 


963 





025 


55123 


225 


16 


9 4 2 





013 


55302 


479 


16 


424 





010 


55303 


411 


16 


247 





012 


55303 


459 


16 


198 





014 


55310 


370 


16 


053 





022 


55310 


413 


16 


164 





010 


55316 


420 


16 


372 





018 


5 5 32 2 


359 


15 


714 





017 


55 32 2 


403 


16 


308 





009 


5 5 328 


453 


15 


988 





Oil 


55335 


477 


16 


124 





010 


55340 


429 


16 


140 





008 


5531 


474 


16 


102 





007 


55345 


470 


16 


182 





016 


55346 


270 


16 


058 





015 


55316 


314 


16 


157 





017 


55351 


271 


16 


196 





(111) 


55356 


309 


16 


165 





013 


5 5 3 5 6 


353 


16 


141 





008 


55361 


311 


16 


156 





009 


55361 


355 


16 


163 





Oil 


55366 


318 


16 


046 





012 


55366 


362 


16 


026 





009 


55371 


384 


15 


963 





Oil 


55376 


377 


16 


1 19 





009 


55376 


420 


16 


1 10 





012 


55380 


467 


15 


931 





009 


55381 


405 


16 


I' "' 





009 


55381 


455 


16 


062 





008 


55386 


397 


15 


940 





007 


55386 


444 


15 


933 





006 


55391 


400 


15 


6 6 2 





013 


55391 


444 


15 


705 





Oil 


55396 


411 


15 


123 





007 


55396 


466 


14 


971 





008 


55401 


304 


14 


843 





012 


55401 


348 


14 


822 





013 


55407 


294 


14 


261 





013 


55407 


342 


14 


263 





Oil 


5 5 410 


2 9 5 


14 


057 





007 


5 5 410 


3 10 


14 


015 





009 


55413 


292 


13 


893 





008 


55413 


336 


13 


886 





006 


55416 


287 


13 


755 





009 


55416 


330 


13 


774 





013 


5 5 1 1 9 


"V 


13 


568 





009 


55419 


326 


13 


567 





Oil 


55426 


161 


<13.409 






55429 


177 


<13.444 






5 5 4 2 9 


220 


<13.420 






55435 


303 


<13.297 






55438 


330 


<13.124 






5 5 13 8 


413 


<13.119 






55441 


365 


<13.231 






55441 


430 


<13.201 






5,> 444 


390 


<i: 


.084 






5 5 444 


434 


<i: 


.041 










<i. 


.017 






55448 


189 


<12.990 






55451 


156 


<13.134 






55454 


265 


<13.104 






55460 


339 


<13.190 






55460 


388 


<13.159 






55463 


341 


<13.137 






55463 


386 


<13.128 






55466 


326 


<13.112 






55477 


191 


<13.173 






55477 


234 


<13.157 






55484 


263 


<13.166 






55485 


119 


<13.260 






55485 


163 


<i: 


.299 






55497 


098 


<i: 


.233 






55497 


141 


<13.232 







